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Abstract 



The inclusive forward jet cross section in deep inelastic e'^p scattering has been measured in 
the region of x-Bjorken, 4.5 ■ 10"'^ to 4.5 ■ 10~^. This measurement is motivated by the 
search for effects of BFKL-hke parton shower evolution. The cross section at hadron level as a 
function of x is compared to cross sections predicted by various Monte Carlo models. An excess 
of forward jet production at small x is observed, which is not reproduced by models based on 
DGLAP parton shower evolution. The Colour Dipole model describes the data reasonably well. 
Predictions of perturbative QCD calculations at the parton level based on BFKL and DGLAP 
parton evolution are discussed in the context of this measurement. 
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1 Introduction 



One of the significant discoveries made at HERA was the steep rise of the proton structure 
function F2{x,Q'^) in the region of small x-Bjorken (x < 10~^) 0]. Various attempts have 
been made to predict the behaviour of F2. In one approach the x dependence of F2 is fitted 
at a fixed value of the scale and then evolved taking into account In terms according to 
the DGLAP evolution equations 0, Different starting scales have been chosen for the 
evolution, down to values as small as ~ 0.3 GeV^ Q. In another approach the leading terms 
in ln(l/x) which appear together with InQ^ terms in the evolution equations are resummed 
to yield the BFKL equation The appeal of the BFKL approach is the fact that it directly 
predicts the behaviour of F2 as a function of x. The CCFM equation p[ interpolates between the 
two approaches and reproduces both the DGLAP and the BFKL behaviour in their respective 
ranges of validity. From the existing F2 data it is not possible to determine unambiguously 
whether the BFKL mechanism plays a significant role in the HERA x-range. 

Several exclusive measurements have been proposed in order to find evidence for BFKL effects 
1^, |1^, |Tl], |T2[. The method proposed by Mueller [0, |I^ starts from the observation that 



the BFKL equation predicts a different ordering of the momenta in the parton cascade (see 
figure |ip. The DGLAP equation predicts strong ordering of the parton transverse momenta /ct,? 
while the BFKL equation relaxes this ordering but imposes strong ordering of the longitudinal 
momenta Xf 

DGLAP: X = Xn< Xn-i < ... < xi, = A;|,„ > k^^^^i > ... > k^^i, 

BFKL: X = Xn <^ Xn-i <^ ... ^ Xi, No ordering in kr- 



Thus the BFKL equation predicts additional contributions to the hadronic final state from par- 
tons with large transverse and longitudinal momenta, i.e. high transverse momentum partons 
going forward in the HERA frameQ. These partons may be resolved experimentally as jets and 
result in an enhancement of the forward jet cross section at small x over NLO calculations and 
parton shower models based on DGLAP evolution. 

The HI collaboration has published results on the forward jet cross section and the transverse 
energy flow in the forward direction and found trends compatible with the expectations of 
BFKL dynamics. 

In this paper we study forward jet production in a region which is expected to be sensitive to 
BFKL parton dynamics {E^ ^ Q"^). The comparison with standard QCD-inspired Monte 
Carlo models is also extended to the range where either E^ <ti or E^ ^ Q^. The mea- 
surement of the forward jet cross section is based on an order-of-magnitude increased statistics 
compared to . Comparisons with different Monte Carlo models are presented and theoretical 
predictions are discussed in the context of this measurement. 



^ In this paper we use the standard ZEUS right-handed coordinate system, in which X = Y^Z = Ois the 
nominal interaction point. The positive Z-axis points in the direction of the proton beam, which is referred to 
as the forward direction. The X-axis is horizontal, pointing towards the centre of HERA. 
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Figure 1: Parton ladder diagram contributing to jet production in DIS. The emitted partons 
in the ladder extend from the bottom of the quark box down to the proton. The transverse 
momenta, kx^i, and the longitudinal momentum fractions, Xi, with respect to the momenta of 
the incoming proton, evolve along the ladder. 



2 Experimental Setup and Trigger 

In 1995 HERA operated with 174 colliding bunches of Ep = 820 GeV protons and = 
27.5 GeV positrons. Additionally 21 unpaired bunches of protons or positrons circulated to 
determine the beam-related background. The integrated luminosity taken in 1995 and used in 
this analysis is 6.36 pb^^. A detailed description of the ZEUS detector can be found in [^. In 
the following the detector components relevant for this analysis are briefly described. 

The central tracking detector (CTD) is surrounded by a superconducting solenoid. Outside the 
solenoid is the uranium calorimeter ||16|, which is divided in three parts: forward (FCAL), barrel 
(BCAL) and rear (RCAL) covering the polar angle region of 2.6° to 176.2°. The calorimeter 
covers 99.9% of the solid angle, with holes of 20x12 cm^ in the centre of the rear and of 
20 X 20 cm^ in the forward calorimeter to accommodate the beam pipe. Each of the calorimeter 
parts is subdivided into towers which are segmented longitudinally into electromagnetic (EMC) 
and hadronic (HAG) sections. These sections are further divided into cells, each of which is 
read out by two photomultipliers. The size of the central FGAL EMG (HAG) cells is 20 x 5 cm^ 
(20 X 20 cm^). From test beam data, energy resolutions of as/E = 0.18/^/E for positrons and 
ge/E = 0.35/v^ for hadrons {E in GeV) were obtained. The effects of uranium noise are 
minimised by discarding cells in the EMG or in the HAG sections if they have energy deposits 
of less than 60 MeV or 110 MeV, respectively. 



A three-level trigger is used to select events online |]I5[. The selection criteria are similar to 
those for the measurement of F2 fl^. For greater than 10 GeV^ the trigger efficiency for 
DIS events is above 97%. About 7 ■ 10® events survive our trigger requirements. 
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3 Event Selection 



Events are selected offline if the Z-component of their primary vertex reconstructed from tracks 
in the CTD lies within ±50 cm around the nominal vertex. A positron candidate with an 
energy E^' above 10 GeV is required. Candidates hitting the RCAL front face within a box of 
X X y = 13 X 8 cm^ around the beam line are rejected because of possible energy loss into 
the beam hole. The fractional energy transfer yjB by the virtual photon in the proton rest 
frame, calculated from the hadronic energy in the calorimeter, is required tohe yjB > 0.1. 
This assures sufficient hadronic energy in the calorimeter to measure the event parameters with 
good accuracy. The i/ei, i.e. the y parameter calculated from the positron energy and angle, is 
required to be below 0.8. This rejects photoproduction events with low energy fake positron 
candidates in the FCAL. The quantity E — Pz = T,iEi[l — cos6i) is required to lie between 35 
and 65 GeV. This requirement also removes photoproduction background. Here Ei and 6i are 
the energies and the polar angles of the calorimeter cells. For DIS events in an ideal detector 
E — Pz is equal to twice the energy of the incoming positron. Events within the x range of 
4.5 ■ 10~^ to 4.5 ■ 10^^ are finally selected for the cross section measurement. 



4 Monte Carlo Simulation 



Monte Carlo event simulation was used to correct the measured distributions for detector 
acceptance and smearing effects. The detector simulation is based on the GEANT 
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program 



and incorporates our understanding of the detector and of the trigger. Events were generated 
with DJANGO 6.24 |g, which interfaces HERACLES 4.5.2 ^ to either the Colour Dipole 



model [H] as implemented in ARIADNE 4.08 or to LEPTO 6.5 H. HERACLES includes 
photon and Z'^ exchanges and first order electroweak radiative corrections. All samples were 
generated with the proton structure function CTEQ4D which describes the measured F2 
structure function |1[]. 

The Colour Dipole model treats gluons emitted from quark-ant iquark pairs as radiation from 
a colour dipole between two partons. This results in partons which are not ordered in their 
transverse momenta /ct- Thus ARIADNE is frequently referred to as "BFKL-like" although it 
does not make explicit use of the BFKL equation. The hadronisation in ARIADNE is based 
on the LUND string model as implemented in JETSET p3 |. 

Another event sample was generated with the MEPS option of DJANGO, as implemented 
into LEPTO. Here the hard interaction is taken from the first order matrix element, but the 
higher orders are simulated by a parton shower based on the DGLAP equation. Thus this 
sample provides a /c-r^ordered parton shower. The hadronisation of the partons is done in the 
same way as in ARIADNE, i.e. via JETSET. In addition to the hard processes, this version 
of LEPTO has non-perturbative effects implemented, so-called Soft Colour Interaction (SCI), 
which affect the investigated phase space region. 



A third sample was generated with HERWIG 5.9 p6|, which like LEPTO has a parton shower 
evolution based on DGLAP. However, there are some differences with respect to LEPTO, such 
as the implementation of colour coherence effects or gluon splitting in the cluster hadronisation 
of the fragmentation phase. 
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A fourth sample was generated at hadron level only with the Linked Dipole Chain model 
option (LDC, version 1.0) of ARIADNE, which uses its own structure functions. The 
parametrisation of "set A" was used, which fits data from HI and ZEUS. In this model the 
parton shower evolution is based on a reformulation |^ of the CCFM approach which 
approximates the BFKL prediction at low x and the DGLAP prediction in the high-x limit. 

All Monte Carlo events were generated with the default settings of the input parameters. 
The first three samples were passed through the full trigger and detector simulation and were 
analysed in the same way as the real data. Additional Monte Carlo samples were generated 
but not passed through the detector simulation in order to study the parton and hadron-level 
properties of the jets. 

The effect of initial and final state QED radiation was studied by generating ARIADNE events 
with and without QED corrections. After applying all event and jet selection cuts the jet cross 
sections as a function of x with and without QED corrections agree within 5%. We conclude 
that QED radiation effects are small and ignore them in the following. 



5 Jet Finding Algorithm and Jet Selection 
5.1 Jet Algorithm 

The analysis is performed with the cone jet algorithm according to the Snowmass convention 



29| . The algorithm is applied to calorimeter cells in the laboratory frame where the cells 
assigned to the scattered positron are excluded. 

The algorithm maximises the transverse energy flow Ex through a cone of radius 



R = y (A?7)2 + (A0)2. Here i? = 1 is used and Arj and A(f) are the differences of pseudorapidi- 
ties and azimuthal angles with respect to the jet direction. The axis of the jet is calculated as 
the transverse energy weighted mean of the pseudorapidity and azimuth of all calorimeter cells 
belonging to the jet. The transverse energy threshold for the seed cells in the algorithm is set 
to 0.5 GeV. Two jets are merged if the overlapping energy exceeds 75% of the total energy of 
the jet with the lower energy. Otherwise two jets are formed and the common cells are assigned 
to the nearest jet. The massless option, where Ex^jet = Pt jet, is used. 



5.2 Jet Selection Criteria 

In order to account for the energy loss of the jets in the inactive material of the detector an 
energy correction procedure is applied both to the total and transverse energies of the jets. The 
correction functions were obtained from the Monte Carlo simulations and are parametrised as 
a function of the total and transverse jet energy. 

After the energy correction, several cuts are applied to select forward jets. The pseudorapidity 
of the jet is restricted to r^jet < 2.6 (equivalent to Oj^t > 8.5°). In this region the jets are well 
reconstructed. The transverse energy of the jets, measured with respect to the direction of the 
incoming proton, is above 5 GeV. The scaled longitudinal momentum xjet = Pz,jet/820 GeV 
is above 0.036. This selects forward jets. The cut 0.5 < E^jet/Q"^ < 2, together with the 
xjet cut, selects the phase space region where BFKL effects are expected. In about 2% of the 
events, two forward jets survive our selection cuts. In these cases the jet with the largest xjet 
is selected. All these cuts restrict the of the selected events to values above 12 GeV^. 
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For events at large values of x, the jet coming from the scattered quark can go sufficiently far 
forward to survive our cuts. In order to reject these events, all found jets are boosted into 
the Breit frame [Q. The boost is calculated from the four-momentum of the virtual photon, 
which is taken as the difference between the incoming and the outgoing four-momenta of the 
positron. Those jets which are in the current region of the Breit frame, i.e. which have negative 
Z-momentum, are rejected. This only affects the two highest x bins where up to 50% of the 
jets are rejected. 

All the selection criteria relevant for the phase space region which defines our cross section 
measurement are listed in table |I|. A total of 2918 events with forward jets survive these cuts. 



E,, > 10 GeV 
y>0.1 
VJet < 2.6 
Erjet > 5 GeV 
Xjet > 0.036 
0.5 < EIjJQ' < 2 
Pz,jet{Breit) > 
4.5 ■ 10"^ < X < 4.5 ■ 10^2 

Table 1: Selected phase space region for the cross section measurement. 



6 Comparison of Data and Monte Carlo Distributions 

Before applying the jet selection criteria, the data and the Monte Carlo predictions are com- 
pared. They agree well in absolute normalisation and in shape. 

After applying the jet selection cuts the uncorrected differential cross sections as a function of 
the event related quantities Q^, E^', yjB and E — Pz are compared in figure 0. ARIADNE 
describes the measured distributions reasonably well while the LEPTO and HERWIG cross 
sections are too small. 

In figure ^ we show uncorrected detector-level cross sections as a function of the jet-related 
quantities Ex jet, xjet, VJet and E^j^^/Q"^. The distributions are compared to the predictions 
of the various Monte Carlo models. All selection criteria are applied except the one for the 
displayed variable. The data in the shaded areas are excluded from the final cross section 
measurement. ARIADNE describes the data in the first three distributions both in shape and 
in absolute value over the entire range. HERWIG and LEPTO underestimate the cross section 
significantly and also disagree in shape. 

The distribution in figure |^d) can be subdivided into three regions. For small E^ j^^/Q'^, i.e. the 
classical DIS regime, all three models agree with the data. Here is large compared to E^ 
and the DGLAP-based Monte Carlo models are expected to describe the data. In the unshaded 
region, which is selected for this analysis, only ARIADNE reproduces the data. HERWIG and 
LEPTO predict much smaller cross sections. In this area we expect significant contributions 
from BFKL-based parton showers. For higher values of E"^ j^^/Q"^ no model agrees with the 
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data. The cross section of ARIADNE is too large, whereas those of LEPTO and HERWIG are 
too smalL In this regime the hard scale is no longer given by the invariant mass squared, Q^, 
of the virtual photon, but by the E"^ of the jets. 



A study of energy flows similar to |^ has been performed in order to investigate whether jets 
in the forward region of the detector still have a pronounced signature and how the beam hole 
and the proton remnant affect the selected jets. 

In figure | we show the transverse energy flow with respect to the forward jet axis averaged over 
all selected events as a function of A0 and Ar], the difference in azimuth and pseudorapidity of 
the cells with respect to the jet direction. The grey bars indicate energy deposits in cells which 
are attributed to the forward jet. The black bars indicate the contributions from those cells 
situated in the towers directly surrounding the forward beam hole. As can be seen, some black 
bars also belong to the jet. The white bars indicate energy deposits which belong neither to the 
jet nor to the cells surrounding the FCAL beam hole. For increasing rij^t the black band-like 
structure in the forward direction of the jet, which we attribute to the proton remnant, becomes 
more and more prominent. For rjjet > 2.6 the selected jets pick up significant contributions 
from the proton remnant in their tails. Studies of the reconstruction accuracy of the angle 
and of the energy of the jets also show a degradation at rjjet values above 2.6 (not shown). 
Therefore, we require rjj^t < 2.6 for this analysis. 

In figure |^ we show the integrated jet shape "^oeti^), i-e. the relative amount of transverse 
energy deposited inside a cone of radius r < R with respect to the jet axis. This function is 
defined as 

where ET{r) is the sum of the transverse energies of all cells of a given jet within a radius r 
with respect to the jet axis. ARIADNE describes the distribution well for all values of Ex^jet- 
LEPTO generates broader jets than observed in the data. The jets are more coUimated as 
Ej-^jet increases. A similar level of agreement between the data and the Monte Carlo events is 
found when bins of rjjet instead of Er^jet are investigated (not shown). 



7 Jet Finding Efficiencies and Purities 

A detailed study of the jet reconstruction quality was performed in order to find acceptable 
cuts for the analysis. 

ARIADNE was used for the study of the efficiencies and purities of the jet reconstruction and 
of the acceptance correction since it describes the data best in shape and absolute rate. The 
efficiency e and purity p of the jet finding are determined as a function of x and are defined by: 

Number of jets'^'^^®^"''^ Number of jets'^^*®^'*'^ 

Number of jets'^"-'^ ' Number of jets''-^^ 

The indices det and had correspond to jets found at the detector or hadron level, respectively. 
Hadron-level jets are defined as jets found by the jet algorithm when applied to the stable 
hadrons from the event generator. The symbol det®had means that the jet has to be found 
at both levels in the relevant variable range and in the same x-bin. The det jets are those 
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Figure 2: Uncorrected detector-level cross sections after the event and jet selection cuts as a 
function of a) Q^, h) E^i, c) y, d) E — Pz- Data are shown as points, ARIADNE as the 
full histogram, LEPTO as the dashed histogram and HERWIG as the dotted histogram. Only 
statistical errors are shown. 
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Figure 3: Uncorrected detector-level cross sections are shown as a function of a) ET^jet, b) Xjet, 
c) VJet, d) E\j^^lQ'^. All event and jet selection cuts are applied (see text) apart from the cut 
on the displayed variable. Data are shown as points, ARIADNE as the full histogram, LEPTO 
as the dashed histogram and HERWIG as the dotted histogram. The data in the shaded areas 
are excluded from the final cross section measurement. Only statistical errors are shown. 
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Figure 4: The transverse energy flow around the forward jet axis averaged over all selected 
forward jets for various rjjet regions. Grey bars indicate the transverse energy in the calorimeter 
cells attributed to the jet, black bars correspond to the transverse energy in the cells around the 
beam hole, and white bars correspond to the transverse energy elsewhere in the calorimeter. 
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Figure 5: Uncorrected integrated jet shapes in intervals of Er^jet which are chosen to contain 
approximately equal numbers of jets. Data are shown as dots. Predictions from ARIADNE 
(LEPTO) are shown as full (dashed) lines. All event and jet selection cuts are applied. The 
statistical errors are negligible compared to the size of the dots. 



surviving all the event and jet selection criteria and the had jets have to survive only those 
cuts which define the phase space region for which the final forward jet cross section is given, 
see table |I|. Figure ^ shows as a function of x the efficiencies e, purities p and correction factors 
c = p/e for the acceptance correction from detector to hadron level. The x bins are chosen such 
that the bin width is at least 2-3 times the x resolution and the statistical errors are below 
20%. The values for e and p lie between 20% and 50%, while those of the correction factors are 
between 1.0 and 1.5. The drop in e and p for x > 10^^ is due to the degraded resolution of x in 
this region of x. The small overall values of e and p are a result of the jet selection cuts and of the 
finite resolutions of the jet variables. The resolutions are: AExjet/ Exjet — 10%, Arjj^t — 0.1, 
/^xjet/xjet ^ 11%, and A{E^ j^JQ^) / {E'^ jjQ'^) ~ 25%. The latter has the largest effect 
on e and p. When the cut on E'^j^^/Q'^ is dropped the efficiencies and purities increase by 
about a factor of two. The measured detector-level rates are multiplied bin-by-bin using the 
correction factor c in order to obtain the hadron-level distributions. In an independent analysis 
the acceptance correction was evaluated using the Bayes unfolding method [^. This method 
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Figure 6: Purity, b) efficiency and c) correction factors for the correction of the detector-level 
jet cross section to the hadron-level jet cross section as a function of x. 

takes migration effects between tlie bins properly into account and is a cross ctieck of tlie 
reliability of the bin-by~bin acceptance correction. The extracted forward jet cross sections 
agree well between the two methods. 



The analysis was also repeated using the clustering algorithm |33| in the Breit frame with 
as scale and a resolution parameter y^ut = 0.5. With these settings the algorithm creates a 
large number of jets at the detector level, which do not have a corresponding jet at the hadron 
or parton level. A change of the ?/cnt^parameter does not improve this situation. Purities and 
efficiencies in the lowest x bins are therefore very small, around 5% and 15%, respectively. 
Nevertheless, the fcr analysis leads to conclusions consistent with those drawn in this paper. 



8 Comparison of DGLAP and BFKL Approaches 

Perturbative QCD predicts the dynamics of the parton evolution. The conventional method 
to solve the parton evolution equations is the DGLAP approach |^, which resums leading 
order (LO) terms proportional to (InQ^)". This results in a parton cascade strongly ordered 
in transverse momentum /c-p, where the parton with the highest transverse momentum appears 
at the lepton vertex of the chain. The longitudinal momenta xj^t decrease towards the photon 
vertex. 

Next-to-leading-order (NLO) calculations, i.e. full second order matrix element calculations 
including first order virtual corrections, where parton densities are incorporated according 
to the DGLAP scheme, are available in three program packages MEPJET, DISENT and 



DISASTER++ pi, RSl BUf. The programs use different techniques to calculate cross sections 



but nevertheless agree reasonably well in their predictions [^|. These NLO calculations are 



not available in full Monte Carlo event simulations. They are purely parton-level calculations, 
delivering parton four momenta which can be analysed, such that the parton level cross sections 
can be evaluated using different jet algorithms and recombination schemes. 
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The BFKL approach of the parton evolution resums terms proportional to (Inl/x)", which 
become dominant over the InQ^ terms at small x. This approach is expected to be valid in 
the high energy limit, where the total available energy, W, is large with respect to any other 
hard scale, Exjet or Q, in DIS. The first term of this resummation is second order in the strong 
coupling constant and is therefore included in the next-to-leading order tree-level diagrams 
in DGLAP-based calculations, e.g. in MEPJET. In figure |l] this term corresponds to exactly 
one parton rung in the gluon ladder between the quark box and the proton. In the following 
it will be referred to as the BFKL 1** term. Since the present approach is only leading Inl/x, 
the parton emissions are strongly ordered in xjet- Recent calculations of next-to-leading Inl/x 
terms in the BFKL kernel |^ predict large negative corrections due to a weakening of the 
strong ordering in xjet which are expected to reduce the cross section significantly. 

The present BFKL calculations do not allow the implementation of a jet algorithm. Therefore 
this calculation can only be regarded as approximate, since the measured jet rates depend on 
the jet algorithm and on their resolution parameters, scales and recombination schemes. For 
example, in our selected phase space region MEPJET yields cross sections for the cone and the 
krp algorithm with their particular choice of parameters which differ by up to 15%. 

In figure |^ we compare the differential forward jet cross section prediction from MEPJET to 
the BFKL 1'^* term and to the leading order (LO) BFKL calculation. We have applied the cone 
algorithm within MEPJET, the same algorithm as used for the data. The renormalisation scale 
in the MEPJET program is varied between 0.25 and 2 to study the scale dependence of 
the parton-level cross section. Here fc^ is the scalar sum of the transverse momenta of the jets 
in the Breit frame. The result changes by ~30% in the two small-x bins and less than 10% in 
the other bins. This is indicated by the shaded band in figure |^. 

The MEPJET NLO and the BFKL l**-term calculations are similar and predict a much smaller 
cross section than the LO BFKL calculation, which shows a steep rise towards smaller values 
of X. 

Also shown are the parton-level cross sections from LEPTO, HERWIG, ARIADNE and from 
LDC Both the MEPS-based LEPTO and HERWIG models show reasonable agreement with 
the MEPJET calculations, whereas ARIADNE exhibits a stronger increase of the cross section 
for small x. The LDC model is well below the ARIADNE predictions. For increasing x all 
models and calculations converge. 

For a direct comparison of the data to the theoretical calculations, the measurements need to be 
evaluated at the parton level, where partons are defined in the Monte Carlo at the stage after the 
last branching of the parton shower and before the hadronisation. The size of the corrections 
from hadron to parton level is studied using the Monte Carlo simulation programs and are 
displayed in figure ^. ARIADNE yields factors close to unity and shows no dependence as a 
function of x. LEPTO and HERWIG show large corrections for small x values. This is expected, 
because these DGLAP based models, which have LO matrix element calculations implemented, 
can only produce a significant number of forward jets due to hadronisation effects and detector 
smearing. The LDC corrections are intermediate to LEPTO/HERWIG and ARIADNE. As 
shown above LEPTO and HERWIG also fail to describe the cross section both in absolute 
size and in shape. Furthermore, the relation between the parton level in parton shower Monte 
Carlo programs and partons in exact NLO calculations is not obvious. Therefore we refrain 
from quoting measurements corrected to the parton level. 
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Figure 7: Parton-level predictions for forward jet cross sections as a function of x. The LO 
BFKL calculation is indicated by the full curve and the BFKL 1^^-term calculation as the 
dashed-dotted curve. The shaded band gives the range of results obtained with MEPJET using 
a renormalisation scale between 0.25 k"^ and 2.0 k"^. The predictions from ARIADNE (full 
histogram), LEPTO (dashed histogram), HERWIG (dotted histogram) and LDC (thin full his- 
togram) are also shown. 




Figure 8: Correction factors corresponding to the ratio of the parton level to the hadron-level 
forward jet cross section as a function of x. The factors for ARIADNE are shown as full circles, 
for LDC as open circles, for HERWIC as open rectangles and for LEPTO as full rectangles. 
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9 Systematic Studies 



We have studied the effects of the variation of several selection cuts and reconstruction un- 
certainties on the final cross section. Figure ^ shows the percentage change of the final cross 
section in each x interval for the major contributions to the overall systematic error. 

Change of the E - Pz cut from > 35 to > 40 GeV 

This tests the amount of photoproduction background in the sample and changes the result by 
less than 6%. 

Alignment uncertainty between the CTD and the FCAL 

The primary event vertex is determined from tracks in the CTD. In order to account for 
an alignment uncertainty between CTD and FCAL, the Z-position of the vertex is shifted 
by ± 0.4 cm. This affects mostly parameters calculated for the scattered positron and related 
quantities like x, and the boost to the Breit frame. The uncertainty from this effect is around 
5%, except in the highest x bin where it reaches 14%. Here, due to the misreconstruction of 
the kinematic variables the current jet may be reconstructed sufficiently far forward to survive 
our selection cuts. 

Jet energy scale uncertainty 

The energy of the jets is scaled by ±5% in the Monte Carlo, reflecting a global uncertainty of 
the hadronic energy scale in the forward region of the FCAL. The result changes by less than 
15%. 

Electromagnetic energy scale uncertainty of the calorimeter 

The energy of the scattered positrons in the RCAL is scaled by ±1% in the Monte Carlo 
corresponding to the global uncertainty of the electromagnetic energy scale in the calorimeter. 
The result changes typically by less than 5%. 

Uncertainty from jet selection criteria 

A possible mismatch between the distributions of jet variables in the data and in the Monte 
Carlo is tested by a variation of the cut values by about one sigma of their resolution followed by 
the determination of the cross sections at the nominal value of the cuts at the hadron level. The 
tested cuts are: the minimum Ex^jet (changed from 5 GeV to 4.5 and 5.5 GeV), the minimum 
xjet (changed from 0.036 to 0.042 and to 0.030), the maximum r]jet (changed from 2.6 to 2.7 
and 2.5), and the minimum and maximum E\IQ^ (changed from 0.5 to 0.6 and 0.4 or from 2 to 
2.4 and 1.6, respectively). All these effects are at the level of a few percent, except in the lowest 
X bin, where they add up to 18%. They are combined (i.e. added in quadrature separately for 
the positive and negative changes) and shown as "jet cuts" in figure ^ 

Acceptance correction with LEPTO 

The full acceptance correction is performed with LEPTO instead of ARIADNE. Since the 
^Tjet/Q"^ distribution of LEPTO differs substantially from the measured one, the LEPTO 
events were reweighted to reproduce the observed E'^ j^^/Q"^ distribution. This test changes the 
cross section by less than 15%, except in the lowest and the highest x bins, where it increases 
the result by +20% and +60%, respectively. 

Further checks included the variation of the accepted ranges of the primary event vertex 
(±40 cm or ±60 cm instead of ±50 cm), a change of the region for the rejected scattered 
positrons within a box of 13 x 8 cm^ or 14 x 14 cm^ around the RCAL beam pipe, and a 
variation of the Uei cut from < 0.8 to < 0.95. The effect on the cross section is negligible. A 
global uncertainty of 1.5% coming from the luminosity measurement is not included. 
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Systematic errors - hadron level 
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Figure 9: Percentage change of the forward jet cross section due to various systematic checks 
in the different x bins. The vertical lines indicate the statistical errors in each interval. Empty 
bins contain values outside the displayed range (see text). 
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10 Hadron— Level Forward Jet Cross Section 



In figure ^ we present the hadron-level forward jet cross section. Tlie numerical values are 
given in table |^. ARIADNE describes the hadron-level forward jet cross section reasonably 
well, apart from the small-x region, where it is slightly below the data. LEPTO and HERWIG, 
as well as LDC, predict significantly smaller cross sections. Switching off the Soft Colour 
Interaction in LEPTO decreases the cross section in the smallest x bin by about 50%, but does 
not affect the large-x region. Increasing the probability of Soft Colour Interaction from the 
default value of 50% to 100% does not increase the cross section. 

The region of large x can be seen more clearly in figure p!0|b). The data, ARIADNE and 
LEPTO converge at larger x. In this region, where Xjet approaches x, the phase space for 
parton emission is small. Therefore, the cross section is expected to be largely independent of 
the parton shower mechanism. On the other hand, HERWIG and LDC stay below the data. 

The excess of forward jets at small x observed in the data with respect to LEPTO and HERWIG 
may be interpreted as an indication of hard physics not implemented in present models of 
DGLAP-based parton evolution. However, the current implementation of BFKL-type physics, 
as exemplified by the LDC model, still underestimates the measured forward jet cross section. 



X range 



da 



II ± Stat. ± syst. [nb] 



syst. {EcAL scale) [nb] 



4.5 ■ 10"^ - 8.0 ■ lO^'* 
8.0 ■ 10-^ - 1.4-10-3 

1.4 ■ 10-3 - 2.5 ■ 10-3 

2.5 ■ 10-3 - 4.5 ■ 10-3 
4.5 • 10-3 - 8.0 ■ 10-3 
8.0 ■ 10-3 - 1.4-10-2 

1.4 ■ 10-2 - 2.5 ■ 10-2 

2.5 ■ 10-2 - 4.5 ■ 10-2 



114 ±9.7 + 29 
96.2±6.5_8.2 



15 
8.2 



77.8±4.7tli 



34A±2.2tii 

14.1 ± i.otii 

6.53±0.54l|]:^ 
2.65±0.25;gi 
0.65 ± 0.09 



(-5.9, +18) 

(-8.1, +7.8) 

(-4.2, +7.0) 

(-2.1, +2.6) 

(-1.2, +1.3) 

(-0.7, +0.2) 

(-0.03, +0.2) 

(-0.00, +0.05) 



Table 2: Forward jet cross sections and their errors for the kinematic region given in table [J. 
The last column shows the correlated systematic error due to the energy scale uncertainty of the 
calorimeter, which is not included in the central column. It corresponds to the shaded band in 
figure [T^. 
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Figure 10: Forward jet cross section at hadron level as a function of x in the kinematic region 
rjjet< 2.6, xjet> 0.036, 0.5 < E^ jJQ^ < 2, Er^jet > 5 GeV, E^, > 10 GeV, y > 0.1. a) Linear 
scale, b) logarithmic scale. Statistical errors are shown as thick error bars, and statistical and 
systematic errors added in quadrature as thin error bars. The errors due to the uncertainty of 
the jet energy scale are shown as the shaded band. 
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11 Summary and Conclusions 



An investigation of forward jet production including a comparison to various Monte Carlo 
models has been performed. Three regions are identified in the j^^/ distribution: i) the 
standard DGLAP region with E^ j^^ <S Q^, where all Monte Carlo models are in agreement 
with the data; ii) the region of phase space where BFKL dynamics is expected to contribute 
significantly with E^ ^ Q^, where only the Colour Dipole model describes the data well, 
and iii) the region with E^ ^ Q^, where none of the models reproduces the data. 

The forward jet cross section at hadron level is measured in the region ii) where E^ ^ 
Q"^. The cross section is compared to the predictions of several models: ARIADNE, which 
includes one of the main features of the BFKL-based phenomenology, that is the absence of 
the strong ordering in the transverse momenta in the parton shower; LDC, which is based on 
the CCFM approach and which smoothly interpolates between the BFKL and the DGLAP 
predictions in their range of validity; and, LEPTO and HERWIG, which are based on leading 
order DGLAP parton evolution. The measured cross section is reasonably well described by 
ARIADNE while LEPTO, HERWIG and LDC predict cross sections that are too low at small 
X. The excess of forward jets at small x observed in the data with respect to LEPTO and 
HERWIG may be interpreted as an indication of hard physics not implemented in present 
models of DGLAP-based parton evolution. However, the current implementation of BFKL- 
type physics, as exemplified by the LDC model, still underestimates the measured forward jet 
cross section. 
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